Spatial nonlocal pair correlations in a repulsive ID Bose gas 
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We analytically calculate the spatial nonlocal pair correlation function for an interacting uniform 
ID Bose gas at finite temperature and propose an experimental method to measure nonlocal corre- 
lations. Our results span six different physical realms, including the weakly and strongly interacting 
regimes. We show explicitly that the characteristic correlation lengths are given by one of four 
length scales: the thermal de Broglie wavelength, the mean interparticle separation, the healing 
length, or the phase coherence length. In all regimes, we identify the profound role of interactions 
and find that under certain conditions the pair correlation may develop a global maximum at a finite 
interparticle separation due to the competition between repulsive interactions and thermal effects. 

PACS numbers: 03.75.Hh, 05.30.Jp, 68.65.-k 
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The study of two-body and higher-order correlations 
is becoming an important theme in the physics of ul- 
tracold quantum gases. Correlation functions are ob- 
servables that provide information about quantum many- 
body wave functions beyond the simple measurement of 
density profiles. They are of particular importance for 
the understanding of low dimensional and strongly cor- 
related systems, atomic gases with exotic phases, and 
systems with multiple order parameters. Such correla- 
tions can manifest themselves in momentum space, as 
with pair correlations of a Fermi gas and can be 
observed in time-of-flight experiments Q which has 
been the standard method of measurement of degener- 
ate gases. Nonlocal correlations in position, on the other 
hand, should ideally be performed in-situ with a spa- 
tial resolution that is smaller than the typical correlation 
length, requiring spatially resolved in-situ single atom 
counting 0, H, 

In this paper we address the problem of nonlocal 
two-particle correlations in a ID Bose gas and propose 
an experimental method to measure them in-situ us- 
ing spatially resolved Raman transitions ^T,] and single- 
atom counting in an optical box trap Q. We calculate 
the spatial second-order correlation function g^'^'^{r) = 
(^t(0)>I'^(r)^(r)*I'(0))/n^ for a uniform gas with repul- 
sive 5-function interactions [i, [IBl (where * is the field 
operator, r is the interparticle separation, and n is the 
ID linear density) and assess the feasibility of the method 
with respect to the characteristic correlation lengths. 

Given a particle at a certain location, the pair correla- 
tion g^^-'(r) describes the probability of finding a second 
particle at a distance r compared to uncorrelated parti- 
cles, and gives the characteristic length scale over which 
the density-density fiuctuations decay. The knowledge 
of g(2)(r) in ID Bose gases is of fundamental impor- 
tance for the understanding of second-order coherence 
and for practical applications such as intensity interfer- 



ometry in ID environments. The recent experimental 
realizations of ultracold atomic gases in the ID regime 

HI, El, [13, [3, El and the fact that the ID Bose gasprob- 
lem is exactly integrable using the Bethe ansatz [3, [lol | 
make this an ideal system for investigating quantum 
many-body physics in previously unattainable regimes. 

An early experimental measurement of atom-atom cor- 
relations using microchannel plate detectors was per- 
formed in an ultracold (but not degenerate) cloud of 
metastable neon [l|. More recently, the method was 
applied to quantum degenerate samples of helium atoms 
and to correlations resulting from condensate collisions 

17| . Other experimental techniques to access higher- 
order correlations include shot-noise spectroscopy of ab- 
sorption images the measurement of three- body 
recombination and photoassociation rates [H, [13, [l9| - 
fluorescence imaging [41, a nd atom counting using high- 
flnesse optical cavities [20[. 

Certain aspects of pair correlations in a repulsive ID 
Bose gas have been studied previously, including the local 
correlation ^^^^(0), asymptotic properties at larg e r and 



25, m 



the zero temperature behavior [21|, [22, l23|, | 
27, bsl bsi. Here we extend these results to nonlocal 



correlations in six analytically tractable regimes of Ref. 



26[ , ranging from strong to weak interactions - all at 
flnite temperatures. 

We begin by recalling that the second quantized Hamil- 
tonian of the system is given by 
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2m 



dx¥¥M, (1) 



where m is the mass and g > is the coupling constant 
that can be expressed via the 3D s-wave scattering length 
a as (7 ~ 2h^a/{ml\) = 2hhj^a [sO]- Here, we have as- 
sumed that the atoms are transversely confined by a tight 
harmonic trap with frequency uj±_ and that a is much 
smaller than the transverse harmonic oscillator length 
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l± = ^/h/nuJ±. The ID regime is realized when the exci- 
tation energy huj± is much larger than the thermal energy 

M. A 



T (with fcs = 1) and the chemical potential [2f 
uniform system in the thermodynamic limit is completely 
characterized 0, [ifl] by two parameters: the dimension- 
less interaction strength 7 = mg/{fi?n) and the reduced 
temperature r — T/Td [25|, where Td = h^ii?/{2m) is 
the temperature of quantum degeneracy. 

Although the uniform ID Bose gas problem is exactly 
solvable by the Bethe ansatz , the cumbersome nature 
of the eigenstates restricts the straightforward calcula- 
tion of correlation functions 22, l32l. The Hellmann- 



Feynman theorem and the solutions to Lieb-Liniger [oil 
or Yang- Yang integral equations - used for calcu- 



lating the local correlation g^'^^ (0) (2^, |25!, |26| - can no 
longer be applied to g^'^\r) at arbitrary separation r. For 
sufficiently large r, the Luttinger liquid theory predicts 
universal features of g^'^\r) [33|. At T = its behavior 
is characterized by interaction-dependent power-law ap- 
proach to the uncorrelated value .g'^^-'(r) = 1 at r ^ cxd. 
At finite T the approach to g'''^\r) ~ 1 becomes expo- 
nential. The Luttinger picture is limited to temperatures 
smaller than the high energy cut-off of the order of Td- 
To describe non-universal features of g^^-* (r) at high tem- 
peratures and short-distances it is necessary to adopt al- 
ternative theoretical techniques [2J, [251, [SJ] . 

Strongly interacting regime [7 ^ max(l, -^/r)]. - We 
employ perturbation theory with respect to the small pa- 
rameter 7~^ [sHi around the Tonks- Girardeau (TG) limit 
of impenetrable (hard-core) bosons 2l| . At T = we ob- 
tain [341 the known |22|. 12911 result 
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where z = nnr. The last term here diverges logarithmi- 
cally with z and can be regarded as a first order pertur- 
bation correction to the fermionic inverse square power 
law. Accordingly, Eq. ^ is vahd for z <C exp(7~^). 

Well below quantum degeneracy, r ^ 1, finite temper- 
ature corrections are obtained using a Sommerfeld expan- 
sion around Fermi-Dirac distribution for quasi-momenta 
at T = 0. For rn <C r^^ this gives [sl] an additional con- 
tribution sin^(7rnr)/127r^ to the rhs of Eq. ([2]), which 
is negligible compared to the T = result as r ^ 1. At 
r = 0, Eq. ^ gives perfect antibunching (7^^^(0) — 0, 
which corresponds to a fully "fermionized" ID Bose gas, 
where the strong interatomic repulsion mimics the Pauli 
exclusion principle for intrinsic fermions. By extending 
the perturbation theory to include terms of order 7"^ we 
reproduce the known T = result for the local correla- 
tion, 5(2) (0) ^ 47rV372 [2i,|25|. 



In Fig. [IK a) we plot the function (/'^'(r), Eq. for 
various 7. The pair correlation exhibits oscillatory struc- 
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FIG. 1: Pair correlation (/'^'(r) as a function of the relative 
distance r (in units of 1/n) in the strongly interacting regime, 
7 3> 1: (a) low temperature TG regime, Eq. ([2]), for r = 0.01; 
(b) regime of high-temperature "fermionization" , Eq. ([Sj. 



ture (similar to Friedel oscillations of a ID electron gas 
with an impurity [si'l), with the local maxima implying 
the existence of more likely separations between parti- 
cles. Despite uniform density, this can be interpreted as 
a quasi-crystalline order (with a period ~ ^/n) in the 
two-particle sector of the many-body wave function. 

Well above quantum degeneracy, r 3> 1, we use per- 
turbation theory combined with the Maxwell-Boltzmann 
distribution of quasi-momenta. In this high-temperature 
"fermionization" regime [25| the characteristic momen- 
tum of the particles is 1/ At and the perturbation param- 
eter is Aid/At <C 1, where aiD = /rag ~ l^ja ^ 
is the ID scattering length and At = ^/^-kW- jmT is the 
thermal de Broglie wavelength. This implies r ^ 7^^, and 
to first order in 7"^ we find js^j 
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(2) 



(r) = 1 - (1 - 2Tnr/7) g-^''-'/^. 



(3) 



In the limit r ^ this leads to perfect antibunch- 
ing, ^'■^•'(O) — 0, while the corrections [as in Ref. [2^, 
(;(2)(0) — 2tI^'^\ are reproduced at second order in 7"-^. 
The correlation length associated with the Gaussian de- 
cay is given by At = \J ^■k j (jn^') . For not very large 7, 
the correlations show a nonmonotonic behavior with a 
global maximum at rmax — 7/2™. This originates from 
the competition between the interaction induced repul- 
sion at short range and thermal bunching [g^^^ (r) > 1] at 
r '-^ At- As 7 is increased the position of the maximum 
diverges and its value approaches 1 in a non-analytical 
fashion g(2)(j.^^^) ~ 1 + (4^/^2) gxp(_^2/8^-)^ 

Figure[T](b) shows a plot of Eq. ([3]) for various 7 and r. 
For a well-pronounced global maximum, moderate values 
of 72/T are required (such as ^"^jr ~ 5, with r = 8, 
7 = 6), and these lie near the boundary of validity of 
our approximations (r^/r 3> 1). The exact numerical 
calculations of Ref. [27[ provide further support for this 
result and show a similar maximum for 72/r ~ 0.25 (with 
r = 47r X 10 and 7 = 10). 

Weakly interacting regime [t^ <C 7 ^ 1] . — For weak 
interactions we rely on the fact that the equilibrium state 
of the gas is a quasi-condensate with suppressed density 
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FIG. 2: Pair correlation g^^\r) in the weakly interacting and 
nearly ideal gas regimes, (a) Solid lines - low-temperature 
weakly interacting gas at r ^ 7 ^ 1, Eq. (j4|; dashed lines - 
DC regime, Eq. ((7}. (b) Solid lines - weakly interacting gas 
at 7 <c: r ^ y/j, Eq. ([5]); dashed lines - DQ regime, Eq. ((U. 



fluctuations and a fluctuating phase [371, l38| . The pair 
correlation function is close to one and the deviations can 



be calculated using Bogohubov theory [2J, |25| . 

At sufficiently low temperatures, r <C 7 <C 1, when 
vacuum fluctuations dominate the excitations and ther- 
mal fluctuations are a smaU correction we find 1341 



■ cosech^ 



2Try^n-^r-^ 87'^/^ 



/ Tmrr \ 



(4) 



where L_i(a;) is the modified Struve function and Ii{x) 
is a Bessel function. The correlation length here is set 
by the healing length ^ — h/^mgn — l/^^n. For r 3> 
^ and finite r, the last term in Eq. ^ dominates the 
others and gives an exponential decay to the uncorrelated 
value of g^'^\r) = 1 (for r one has a power law 
decay). Even at T = 0, oscillating terms are absent, in 
contrast to the strongly interacting regime, Eq. ([2]). The 
limit r — > reproduces the result of Eq. (9) of Ref. [isl . 
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(2) 



(0) = 1 - + 7rTV(2473/2). In Fig. [Ifa) we plot 

Eq. for different values of 7, and we note that the 
finite temperature correction term is negligible here. 

In the opposite limit, dominated by thermal fluctua- 
tions corresponding to 7 ^ t <C -^7, we find [l^] 



?(2)(r) = 1 + 



2V7 



-V7[L-i(2V7^^)-/i(2V7"r)], (5) 

valid for r/^ < 1. The last two terms are due to vacuum 
fluctuations and are a negligible correction, so the leading 
term gives an exponential decay [see Fig. (Hb)] with a 
correlation length given by the healing length ^. The 
peak value is g^^^O) = 1 + T/{2y/j), in agreement with 
Ref. [2^. For r/f ^ 1 vacuum fluctuations dominate 
and we reproduce the asymptotic behavior of Eq. (j4|) . 

Nearly ideal gas regime [7 <C minjr^, -x/t}]- — Finally, 
we present the results for the decoherent regime, where 
both the density and phase fluctuations are large and the 



local pair correlation is close to the result for noninteract- 
ing bosons, .g'-^-'(O) = 2. Depending on the temperature 
r, we further distinguish two sub-regimes: decoherent 
quantum (DQ) regime for r ^ 1, and decoherent classi- 
cal (DC) regime for r ^ 1. Both can be treated using 
perturbation theory with respect to the coupling constant 
g around the ideal Bose gas result of Ref. 3] ■ 

In the DQ regime, with ^/^ ^ r <C 1 [1^ , the nonlocal 
pair correlation is [3J| 



(r) ^ 1 + [1 _ 4^(1 + rnr) /r^ 



(6) 



with the peak value g^^^O) = 2 - A-^/t"^ [251. For 7 = 
the correlations decay exponentially with the character- 
istic correlation length which coincides with the phase 
coherence length ~ h^n/mT = 2/rn [2^ and is re- 
sponsible for the long-wavelength phase fluctuations. For 
7 > 0, g'''^\r) becomes nonmonotonic with a minimum 
at nrniin — t/47 ^ 1 before reaching its uncorrelated 
— > 00) = 1. Thus, at intermediate range 



(2)^ 



value g 

we have weak antibunching due to interatomic repulsion, 
while at short range we have typical ideal-gas bunching 
due to exchange interaction, as shown in Fig. [2Kb). 

In the DC regime (r ^ max{l,72}), the pair correla- 
tion is given by [34] 



1 



7\/27r/r erfc(\/™2r2/2), (7) 



^^'(0) -f 272/T at nonzero separation 
As 7 is increased, there is a con- 



where erfc(x) is the complimentary error function. At 
r = we have g^^HO) = 2 - 7^27r/T [ll] as erfc(O) = 1. 
In the noninteracting limit (7 — 0) we recover the 
well-known result for the classical ideal gas [3§| char- 
acterized by Gaussian decay with a correlation length 
At. For 7 > we observe [see Fig. dKa)] the emer- 
gence of nonmonotonic behavior, with a global maxi 
mum g(2)(,,^^^) ^ g 

nrn^a,^ = 27/T < 1 
tinuous transition from the DC regime to the regime of 
high-temperature "fermionization" , with g^^) (Q'j reducing 
further and the maximum moving to larger distances. 

We now discuss experimental methods to measure pair 
correlations. Local correlations have been measure d by 
photoassociation [l^ and by three-body loss 12, Il9| . 
however, nonlocal spatial correlations have not been mea- 
sured in-situ to the best of our knowledge. We discuss 
one possible implementation of spatially resolved imag- 
ing in the context of a ID gas. The method is closely 
related to the spatially resolved measurement of Ref. [tI] . 
The first step is to "freeze-in" the correlations by turn- 
ing on a deep standing wave of far-detuned light, which 
breaks up the distribution into discrete packets, sepa- 
rated by a half-wavelength. The next step is to apply a 
magnetic field that varies linearly in magnitude, creating 
a spatially dependent Zeeman shift. Next, a stimulated 
Raman transition is driven such that atoms in two sites 
are transferred to a different hyperfine state. This can 
be accomplished by using two frequencies on the Raman 
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beam corresponding to two different locations. The un- 
affected atoms can be removed with a pulse of resonant 
light and the remaining atoms counted. 

For a specific example, consider a degenerate gas of 
sodium atoms in an optical box trap [8|, in the F — 
2,mp = —2 state. After "freeze-in" with a A = 532 nm 
lattice, we apply a magnetic field gradient of 150 G/cm. 
We then drive a two-photon stimulated Raman transition 
to the F — l,mp — —1 state with Raman beams at 532 
nm with optical power of 2.5 mW focused to 50 /im and 
a pulse duration of 260 /zs. A detuning of 53 kHz corre- 
sponds to a shift of one lattice site and therefore a reso- 
lution of ~ 266 nm. For typical ID gas parameters, this 
resolution is sufficient to resolve the characteristic cor- 
relation lengths found here. A related simpler method, 
but with somewhat lower resolution, has been recently 
proposed in Ref. [40| . 

In summary, we have calculated nonlocal pair correla- 
tions in a uniform ID Bose gas in six physically relevant 
regimes. The correlations can be measured using spa- 
tially resolved single-atom counting. We have shown ex- 
plicitly that the characteristic global correlation lengths 
are given by one of four length scales: the thermal de 
Broglie wavelength, the mean interparticle separation, 
the healing length, or the phase coherence length. In all 
cases we identified the profound role of interactions that 
can lead to non-trivial structures with local maxima or 
minima at a finite interparticle separation. 
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